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ABSTRACT 

We investigate the GeV emission from gamma-ray bursts (GRBs), using the results from the Energetic 
Gamma Ray Experimental Telescope (EGRET), and in view of the Gamma-ray Large Area Space Telescope 
(GLAST). Assuming that the conventional prompt and afterglow photons originate from synchrotron radiation, 
we compare an accompanying inverse-Compton component with EGRET measurements and upper limits on 
GeV fluence, taking Klein-Nishina feedback into account. We find that EGRET constraints are consistent with 
the theoretical framework of the synchrotron self-Compton model for both prompt and afterglow phases, and 
discuss constraints on microphysical parameters in both phases. Based on the inverse-Compton model and 
using EGRET results, we predict that GLAST would detect GRBs with GeV photons at a rate > 20 yr _1 from 
each of the prompt and afterglow phases. This rate applies to the high-energy tail of the prompt synchrotron 
emission and to the inverse-Compton component of the afterglow. Theory predicts that in a large fraction of 
the cases where synchrotron GeV prompt emission would be detected by GLAST, inverse-Compton photons 
should be detected as well at high energies (> 10 GeV). Therefore GLAST will enable a more precise test 
of the high-energy emission mechanism. Finally, we show that the contribution of GRBs to the flux of the 
extragalactic gamma-ray background measured with EGRET is at least 0.01% and likely around 0.1%. 
Subject headings: gamma-rays: bursts — radiation mechanisms: non-thermal 



1. INTRODUCTION 

Cosmological gamma-ray bursts (GRBs) have released a 
tremendous amount of energy in the past and present Uni- 
verse. Their emission covers very wide range of frequencies: 
a highly variable prompt phase radiates — 100 keV gamma 
rays, while a subsequent afterglow radiates radio to X-ray 
photons. It is likely that the bulk of these photons are emitted 
by gyration of relativistic electrons in magnetic fields — e.g., 
synchrotron radiation. The relativistic electrons are acceler- 
ated in either internal dissipation (for prompt emission) or ex- 
ternal shocks (fo r afterglows). For reviews, see, Piran (2005); 
lMeszarosl(l2l^ : lN^ail(l2007h . 

GeV photons were detected as well from several GRBs 
by the Energetic Gamma Ray Experimental Telescope 
(EGRET) on board the Compton Gamma Ray Obser- 
yatory (CGRO) ([ Schneid et al. 1992; Sommeretal. 1994; 
iHurlev et alJl994lSchneid et al.ll9 95: Gonzalez et al. 2003). 
The data are still not sufficient for us to firmly infer 
emission mechanisms of these GeV gamma rays, but the 
most promising mechan ism is synchrotron self-Compton 
(SSC ) scattering (e.g ., Meszaro s. Rees. & Papathanassioul 



according to theory (e.g., Sari & Esin 20011) . 

In this paper, we explore the GeV gamma-ray emission 
of GRBs in the context of SSC mechanism. 1 Besides the 
several GRBs detected by EGRET, there are many oth- 
ers for which upper bou nds on the fluence were obtained 
dGonzalez Sanchez! 12005). These -100 GRBs should also 
be compared with the predictions of SSC model, because the 
fluence upper limits in the EGRET energy band are compa- 
rable to the fluence of prompt emission collected by Burst 
And Transient Source Experiment (BATSE) instrument on- 
board CGRO. As the experimental bound is already strong, 
while theoretical models of SSC process predict a large flu- 
ence for the EGRET energy range, we derive meaningful 
constraints from EGRET data analysis on the physics of the 
high-energy emission mechanisms of GRBs. This approach 
is different from (and therefore complementary wit h) that 
in previous stu d ies ( e.g., Dermer, Chiang, & Mitman 2000; 
Asano & Inoue 2007; I okaetalJ l2007t (Gou& Meszaros 
2007; Fan et al. 2008; Murase & Ioka 2008; Panaitescu 2008 
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19941: IWaxmanl Il997t IWei & Lull 19981; ICh iang & Dermer 
Panaitescu & Kumar 2000; Zhang & Meszaros 2001; 



Sari& Esin 2001; Guetta & Granot 2003). This is because 



the relevant emission parameters such as the energy fraction 
of the GRB jets going to electrons (e e ) and magnetic fields 
(eg) are relatively well measured from the afterglow spectra 
as well as light curves; the typical values are e e = 0. 1 an d 
e B = 0.01 (e.g., iPanaitescu & Kumaill200ll lYost et alJl2003h . 
In the prompt emission, e e is similar or even higher, as evident 
from the high efficiency of this phase, while e# is not well con- 
strained. Thus, there should be a significant inverse-Compton 
(IC) component accompanying the synchrotron radiation in 
both the afterglow and prompt emission. The luminosities of 
the synchrotron and IC are expected to be comparable as IC- 
to-synchrotron luminosity ratio is roughly given by (e^/e^) 1 ' 2 , 



and references therein), where the prediction of gamma-ray 
flux relies only on theoretical models and sub-GeV observa- 
tions. We instead use EGRET data in order to infer the GeV 
emission and constrain the theoretical models. 

We use our results to predict the expected number of GRBs 
that would be detected by the Gamma-ray Large Area Space 
Telescope (GLAST). The GLAST satellite is equipped with 
the Large Area Telescope (LAT), which is an upgraded ver- 
sion of EGRET. Since revealing the high-energy emission 
mechanisms of GRBs are one of the important objectives of 
GLAST, our prediction should give a useful guideline. Finally, 
we apply our results to estimate the contribution of GRBs 
to the diffuse extragalactic gamma-ray background (EGB), 

1 Our analysis and conclusions are applicable also if the MeV and/or radio- 
X-ray afterglow emission mechanism is not synchrotron but another type of 
emission from relativistic electrons that gyrate in a magnetic field, such as 
jitter radiation iMedvedev 2000). 
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which was also me asured by EGRET ( Sreekumar et al l 19981: 
Strong et al.1 12004 see, however, [Keshet, Waxman, & Loebl 
2004 for a subtle issue of Galactic foreground subtraction). 

This paper is organized as follows. In § [2] we summarize 
the predict ions of SSC model for the prompt (§ 12. It and af- 
terglow (§ I2.21 i phases. Section [3] is devoted for analysis of 
the GRB fluence data by EGRET, from which distributions of 
fluence in the GeV band are derived. We then use these dis- 
tributions to argue prospects for GRB detection with GLAST 
in §|4] and implications for EGB from GRB emissions in §[5j 
In §[6] we give a summary of the present paper. 

2. INVERSE-COMPTON MODEL OF HIGH-ENERGY EMISSION 

If the prompt and/or afterglow emission is due to syn- 
chrotron radiation from relativistic electrons (with Lorentz 
factor 7 e ), then there must be an accompanying IC compo- 
nent from the same electrons scattering off the synchrotron 
photons. The spectral shape of the IC emission is almost the 
same as the synchrotron radiation (shifted by 7 2 ), and is ex- 
pected to fall around the GeV range during both the prompt 
and afterglow phases. For e e > eb, and assuming that there 
is no "Klein-Nishina suppression" and that the emitting elec- 
trons are fast cooling, the IC fluence is related to the syn- 
chrotron fluence simply through Fie sw (e e /es) 1/ ' 2 F S y n . Thus, 
assuming that the microphysics do not vary much from burst 
to burst, it is natural to assume proportionality between the 
synchrotron MeV fluence (observed by BATSE) and the GeV 
synchrotron plus IC fluence (observed by EGRET and in the 
future by GLAST): 



FG e v = (r] S y n + r} I c)Fu 



(1) 



where ?7 syn and 771c are coefficients for the proportionality due 
to synchrotron and IC processes. Note that the synchrotron 
fluence in the GeV range can be extrapolated relatively eas- 
ily, if we assume that the spectrum extends up to such high 
energies. Thus, we here focus on theoretical evaluation of 
the IC component. At first approximation, the coefficient 771c 
is roughly (e e /e/j) 1//2 from considerations above, and thus we 
define 



, \i/2 F 

?7IC= ( — CKN?w^ 1 



(2) 



where for the prompt emission F syn « FivieV while for the af- 
terglow F syn is the afterglow fluence within the radio to X-ray 
energy bands. Correction factors £kn an d £w represent the 
effect of Klein-Nishina suppression and detector energy win- 
dow, respectively, which are given below. 

We define typical frequencies for both synchrotron (f syn ) 
and IC (t'lc) as the frequencies where most of the energies are 
radiated in case that the Klein-Nishina cross section does not 
play an important role; i.e., where vf v for each component is 
peaked in this case. From relativistic kinematics, these two 
typical frequencies are related through 



V\C ~ 7 m ^syn, 



(3) 



where 7,, is a characteristic Lore ntz factor of the electrons tha t 
dominate the synchrotron power (Rybicki & Lightman 1979); 
this is true in the fast cooling regime, which is the case in 
the most of our discussions (Sari & Esin 2001). The Klein- 
Nishina effect is relevant if a photon energy in the electron 
rest frame exceeds the electron rest mass energy, and this con- 
dition is formulated as 



where Ti, is the bulk Lorentz factor of the ejecta, which is 
on the order of 100 in the prompt phase of GRBs and their 
early afterglows. Upscattering synchrotron photons to ener- 
gies above /wkn is highly suppressed, which results in IC cut- 
off at i/KN- 

Besides producing a spectral cutoff, the Klein-Nishina ef- 
fect also modifies the way electrons cool, which is relevant 
for the GeV emission and is also included in £kn- Elec- 
trons with energies above Klein-Nishina threshold (for a given 
seed-photon energy) can lose their energies only through 
synchrotron radiation, while the lower-energy ones can cool 
through both processes. Such an effect has been studied in the 
case where the seed phot ons for IC scatte ring are provided by 
an external sources (e.g., Moderski et al. 2005a b, and refer- 
ences therein). However, in the case of SSC mechanism, since 
the seed photons are emitted from synchrotron process due to 
the same electron population, we should properly take into ac- 
count feedback. Giving full details on this is beyond the scope 
of the present paper, but some results are sum marized briefly 
in Appendix lAl (see also Derishev et al.ll2003l) . Here we only 
show the approximate analytic form of £kn: 

C 1 for 7,, < 7kn, 

&nH ( 7,„V 1/2 f ^ ( 5 ) 

I (ij for 7« > 7KN, 

where 7kn is the Lorentz factor of electrons for which photons 
at v > i/ syn are in the Klein-Nishina regime. The energy of 
an observed photon with frequency v as measured in the rest 
frame of an electron with Lorentz factor 7 is ss jhu/Ti, where 
the 1/F& factor converts the photon energy from the observer 
frame to the plasma rest frame and the 7 factor converts it to 
the electron rest frame. Since such a photon is in the Klein- 
Nishina regime of an electron with Lorentz factor 7 once its 
energy in the electron rest frame is larger than m e c 2 we obtain: 

F b m e c 2 

7kn = — ■ (6) 

hv syn 

This Klein-Nishina feedback effect modifies the spectrum 
shape of both synchrotron and IC emissions (in addition to 
the Klein-Nishina cutoff for IC). We note that equation <|3J 
provides a solution that agrees within a factor of ~2 with the 
one obtained by numerically solving equation ( lAlb . This pre- 
cision is sufficiently good for our purpose, especially because 
it is well within the uncertainty ranges of other parameters. 

By £i t , we take into account the fraction of the IC fluence 
that falls into the GeV detector energy bands. EGRET win- 
dow is between hv w [ = 30 MeV and hv wu = 30 GeV while 
GLAST-LAlT window is between hv K j = 20 MeV and hu WtU = 
300 GeV. We here assume that the frequency where most of 
the IC energy is released, fic.peak = min[z/ic,J / KN], is always 
larger than lower limit of the frequency band, v w .\, as expected 
for both EGRET and GLAST, and thus consider the cases in 
which z^ic.peak is within or above the detector frequency band. 
In the former case where v w j < i^cpeak < v wu> we h ave £ w ~ 1. 
On the other hand, if ^ic.peak > v KM , then most of the en- 
ergy comes from the upper frequency limit v WM , and we have 
Cw ~ (^H'.u/^iCpeak) 2 " 1 ! where a\ is the photon spectral index 
below peak frequency. Thus we may approximate £ w as 

( , , , mint^cfKNr " : ~" 



hum = T h j m m e c 



(4) 



l + ^Y' ' for 7m < 7KN , 
r "" L,-2 (7) 

1 + T?7-) fo r 7m>7KN, 
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where z^c/^kn = 7»i/7kn as one can easily show. 

The discussion above assumes that the density of the syn- 
chrotron photon field is proportional to the instantaneous syn- 
chrotron emissivity. In the case of relativistically expanding 
radiation front, this assumption is valid when the duration 
over which the emissivity vary significantly, St, is compa- 
rable to the time that passed since the expanding shell was 
ejected, fo. In this case the ratio between the synchrotron 
emissivity and the synchrotron photon field density is in a 
steady state. When St <C fo the synchrotron photon field den- 
sity may be significantly lower than in the steady state case 
dGranot. Cohen-Tanugi. & Siival 120071) . thereby suppressing 
the IC component. The exact suppression factor depends on 
the detailed spatial and temporal history of the emissivity. 
Theoretically, in the afterglow phase we expect St ~ fo. Also 
in the prompt e mission pha se, internal shock models generally 
predict St ~ fo (Pirarj [l999l and references therein). Thus, in 
the internal-external shock model corrections to the IC com- 
ponent due to this effect are expected to be on the order of 
unity. Therefore, in the present paper, we assume that such 
an effect can be neglected and that the synchrotron photon 
field is proportional to the instantaneous synchrotron emis- 
sivity. One should keep in mind, however, that St <C fo is a 
viable possibility (see, e.g.. lPe'er& Waxman 20041 120051 for 
a more detailed study in such cases), especially in the highly 
variable prompt phase. In principle, detailed GLAST observa- 
tions of an IC emission may be able to constrain St /to during 
the prompt phase. 

In addition, towards the higher end of the EGRET 
or GLAST energy band, photons may start to be 
subject to absorption due to pair creation in the 
sourc e or during prop a gation ( e.g.. iBaring & Hardind 
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iLithwick & Saril 1200 ll iRazzaaue et all 
l2004t ICasanoya. Dingus. & Zhangl 
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Although such 



mechanism might be relevant for the IC yields (especially in 
the prompt phase) depending on some parameters that are 
not well constrained yet, we assume that it is not the case in 
the present paper. GLAST will hopefully provide information 
that enables better handle on this issue. 

2.1. Prompt phase 

BATSE (as well as Swift satellite) detected so far a large 
number of GRBs in prompt phase with gamma rays in the en- 
ergy band of 20 keV-1 MeV. The spectrum is well described 
by a smoothly broken power law with a typical lower-energy 
index of a\ ~ 1 and higher-energy index of «2 ~ 2.3; the 
spectral break typically occurs around hv syn « 300 keV, where 
the energy of th e promp t emission vf v pe aks (Ba nd et all 
1993; Pre eceetalJ 12000: Ka neko et al.ll2006l) . As we show 
in Figure Q] the distribution of the fluence integrated over the 
BATSE energy band follows log-normal function. 2 The peak 
of this distribution is Fbatse = 2.5 x 10~ 6 erg cm" 2 , and its 
standard deviation is u\ 0% f = 0.75. The average of the BATSE 
fluence is therefore (Fbatse) = 10~ 5 erg cm -2 . 

Therefore, for the prompt emission phase, using a,\ = 1, 
«2 = 2.3, and hv syn = 300 keV, we find: 

7kn = 170r^ 2 - 



(8) 



where T/,2 = T/,/10 2 . In addition, for £„,, considering GLAST- 
LAT energy window (20 MeV-300 GeV) in equation (|7), we 

2 http://www.batse.msfc.nasa.gov/batse/grb/ 
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FIG. 1 . — The fluence distribution in the prompt phase of BATSE GRBs. 
The best-fit log-normal function is also shown, where the peak is at Fbatse = 
2.5 X 10 -6 erg cirT 2 and the standard deviation is <ri gF - 0.75. 



obtain 



sw : 



i+(»r 



i+ 



5900 



for 7m < 170r 4 , 2 , 
-i (9) 

for 7 „, > 170I\ 2 . 



Now assuming that all electrons are accelerated in the shocks, 
the typical value for the Lorentz factor of the relativistic elec- 
trons are given as 



7» 



e e ^(r re i-l) = 200e e ,_ 1 (r rel - 
m P 



1). 



(10) 



where r re i is the relative Lorentz factor of the colliding ejecta 
portions and t e -\ = e e /10 _1 . In the internal shock model for 
the prompt emission, r, e i - 1 is of order unity. If we adopt 
r re i = 3 and e e = 0.1, we obtain j m rs 400. Furthermore, 
assuming T/, = 100, equation (0 gives £ w « 0.9, and equa- 
tion (0 with equation © gives £kn ~ 0.7. By substituting 
these values and assuming eg = 0.01 in equation (0, we ob- 
tain rjic ~ 1.9, which implies that under the most straightfor- 
ward assumptions a comparable fluence is expected in both 
GLAST-LAlT and BATSE windows. In this case, the Klein- 
Nishina cutoff energy is in GLAST-LAlT band as well as in 
EGRET band (/wkn ^5 30 GeV), and thus we also obtain an- 
other comparable value of 771c ~ 1 .2 in EGRET case. 

Note that in the case of prompt emission, the synchrotron 
spectrum is not negligible in EGRET and GLA5T-LAT en- 
ergy bands. For canonical parameters (hv %ya = 300 keV, 
ot\ = 1, and ct2 = 2.3), the ratio f/„,ic/ l/ ./i/,syn at 100 MeV 
is about 0.01(e e /eB) 1 / 2 (7 m /400)~ 2 , assuming that the syn- 
chrotron spectrum continues into the GeV window without a 
break and IC is not much suppressed by the Klein-Nishina ef- 
fect. Therefore, the synchrotron component dominates around 
the lower-energy limit where most of the photons (although 
not most of the fluence) are observed. In the case of EGRET, 
since only a handful of photons were detected in all EGRET 
events, these are expected to be dominated by the synchrotron 
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low-energy (—100 MeV) photons. This indicates that the 
quantity we can constrain using the EGRET fluence upper 
limits is not 771c but ?7 syn = F syn (100 MeV)/FM e v> the ratio of 
synchrotron fluence around 100 MeV and that in the MeV 
range. In addition, this picture is indeed consistent with the 
fact that the spectral indices of GeV phot ons for several GRBs 
measured with EGRET are a = 2-3 (e.g . , ISchneid et al.ll 19921 
ISommeret"aIlll994t lHurlevetai1ll994l) . Note however, that 
the energy fluence in GLAST- LAT and EGRET bands can be 
dominated by a much harder IC component (a w 1-2) that 
peaks above — 1 GeV and may carry up to — 10 times more 
energy than the one observed at 100 MeV without being de- 
tected. This is because even when the — 10 GeV fluence is ten 
times larger, the small photon number at such high-energies 
is still small enough to avoid detection. Thus, EGRET obser- 
vations, which are consistent with measurement of the syn- 
chrotron high energy tail, can only put an upper limit on rjic- 

2.2. Afterglow phase 

The afterglow is considered to be a synchrotron emission 
from electrons accelerated in the external shock, which is 
caused by the interaction between the relativistic ejecta and 
the interstellar medium. In this model, the synchrotron emis- 
sion dominates the spectrum from radio to X-ray. The asso- 
ciated IC emission is expected to dominate the GeV energy 
range (i.e., 771c 3> f] iyn ), since the electron Lorentz factor is 
much larger than the case of prompt emission (see eq. fTOl , 
where the relative and bulk Lorentz factors are the same, 
T re i = 1$), compensating the smaller v syD (eq. (3))- During 
the first several minutes (observer time), electrons might be 
cooling fast (a\ = 1.5) with hv S yn — 1 keV, while j m sw 10 4 - 
10 5 . This implies that the fraction of the IC energy that 
falls in GLAST-LAlT energy window is close to unity, i.e, 
£ w w 0.2-0.9 from equation © (for EGRET 77,,, « 0.08-0.5) 
and £kn ~ 0.7-1 from equations (O-©. Since hv\c at early 
time is close to the upper limit of the energy window the ef- 
fective photon index of the IC emission within the detector 
window during this time is « 1.5-2. 

At later times the electrons are at the slow-cooling regime 
and r/ syn is the cooling frequency, while a typical j e is the 
Lorentz factor of electrons that cooled significantly (e.g., 
Sari & Esin 2001). In this regime the SSC peak is very broad 
and its location is almost constant with time. For typical pa- 
rameters, the Klein-Nishina effect do not play a major role 
while the peak of the SSC emission falls within GZAST-LAT 
and EGRET windows. Therefore, at late time £ w ss 1 and the 
effective photon index within the energy windows of these de- 
tectors is ?s2. 

One should, however, note that on long time scales the GeV 
background becomes important, making it hard to detect the 
GeV afterglow. Therefore, the optimal time scale for GeV 
afterg low search would be — 100-10 3 s (Zhang & Meszaios 
2001). The afterglow GeV fluence, Fo e v in equation ([T}, is 
that integrated over a given time scale, while Fjviev is col- 
lected over roughly Z90, during which 90% of the MeV pho- 
tons are counted. The total energy radiated away by the ra- 
dio to X-ray afterglow during every decade of time is roughly 
0.01-0.1 of the energy emitted in the prompt phase. There- 
fore we expect a bright GeV afterglow which radiate about 
0.01-0. l(e c /eB) 1 ^ 2 fMeV every decade of time for hours and 
days after the bursts. In this paper when considering EGRET 
observations, we adopt 200 s after 790, when electrons are in 
the fast cooling regime, as the duration over which 7*bev is 
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FIG. 2. — Contour plot of allowed region in /i—tr space, obtained with the 
analysis of EGRET data assuming a spectral index of —2.4 during TgQ (top) 
and 200 s after TgQ (bottom), fi and a are the central value and standard 
devi ation , respectively, for the log-normal distribution of the fluence ratio r\ 
(eq. Hil l. The best fit points (A) are marked as crosses, and other represen- 
tative points (B and C) are also indicated in both panels. 



integrated. 

3. CONSTRAINT ON HIGH-ENERGY EMISSION WITH EGRET 



Gonzalez Sanchez (2005) analyzed GRBs that were de- 
tected by BATSE and observed by EGRET. Since the field of 
view of EGRET was much smaller than that of BATSE and the 
observation was limited by the life time of the spark chamber, 
EGRET covered only about 100 GRBs out of -3000 BATSE 
bursts. But this is still a reasonably large number to get statis- 
tically meaningful result. The analysis of the prompt burst 
in EGRET data was performed around the error circles of 
BATSE bursts for the first 790, and spectral index of -2.4 is 
assumed within EGRET window (the upper limits are higher 
by a factor of «10 for a spectral index of — 1). The same anal- 
ysis was performed f or the afterglow phase, fo r 200 s after 790 
(not including Tgo). Gonzalez Sanchez (2005) measured the 
fluence of 6 and 12 GRBs, in prompt and afterglow phases re- 
spectively. For all other GRBs only fluence upper limits were 
obtained in the range 10~ 6 -10~ 3 erg cm" 2 . 

Here we interpret these results in the framework of the 
SSC model, which implies that the fluences in BATSE and 
EGRET bands are likely to be positively correlated through 
equation (Q) (F B atse = F MeV and F EG ret = T^Gev)- We further 
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FIG. 3. — The efficiency of EGRET for GRBs as a function of fluence, e(F), 
during (a) Tgo and (b) 200 s after Tgo (assuming a spectral index of —2.4). The 
histogram represents cumulative fraction of GRBs whose fluence limits are 
below a given value, which can be interpreted as the detector efficiency, while 
the solid curves are fitting function. 



assume that the coefficient 77 (?7 syn for prompt and rjic for af- 
terglow phases) follows some probability distribution func- 
tion 77(77) which is independent of Fbatse- We consider a log- 
normal distribution with the central value p and standard de- 
viation a: 



p{r]\n,a)dr} = 



1 



27TCT 



exp 



(log 77 -/j) 2 



2a 2 



dlogrj. (11) 



Constraining [i and a then leads to implications of GRB pa- 
rameters such as e e , eg, and j m , through their relations given 
in the previous section. 

We used the observations to constrain [i and a by carry- 
ing out a maximum likelihood analysis. 3 Figure [2] shows the 
contour plot of the most likely region on the p-a plane for 
790 (top) and 200 s after Tgo data (bottom) assuming a spec- 
tral index of -2.4 (if the spectral index is -1 then p increases 
by wl). In that procedure, detection efficiency of EGRET 
as a function of fluence, e(F), is obtained from the distribu- 
tion of the EGRET upper limits (for undetected GRBs), which 
is shown in Figure |3j i.e., a cumulative fraction of bursts 
whose fluence limits are below a given fluence. In the case 
of detected GRBs, on the other hand, the size of the error 
bars for the fluence is interpreted as measurement accuracy of 
EGRET. Then, in order to test the consistency of the assump- 
tion that equation (fTTT i fits the data, we carried out a Monte 
Carlo simulation that draws 10 5 realizations of EGRET ob- 
servations assuming that the distribution of ^egret/^batse 
follows equation ( fTTT i with the most likely values of p and 
a. By comparing the likelihood of these Monte Carlo real- 
izations with that of the actual EGRET observations, we find 
that 70% of the realizations have a lower likelihood, suggest- 

3 The log likelihood of a distribution is calculated by integrating the prob- 
ability between the error bars and below the upper limits of EGRET ob- 
servations. For the Tgp fluence data, we used the results of Fig. 2.3 of 
Gonzalez Sanchez (2005) rather than Tables 2.1 and 2.2 there. 



FIG. 4. — Distribution of EGRET fluences during (a) Tgo and (b) 200 s after 
Tgg. Models A-C correspond to the points on fjr- a plots in Fig. [5] The BATSE 
fluence distribution is also plotted for comparison. The distribution for the 
prompt phase (a) is for the high-energy tail of the synchrotron radiation. The 
prompt IC fluence may be larger by up to about one order of magnitude (see 
text). 



ing that equation ( fTTT i with its most likely values is indeed 
consistent with the observations. 

Given p and a, we can obtain the distribution of fluence 
in EGRET band by convolving BATSE fluence distribution 
(dN/dFsATSE, Fig.fJJ and p(f]\p, a): 



— =/ dr\ p(ri\n,<r) — 

"Regret Jo «-tbatse 



(12) 



17 l Fw 



As representative models, we use three sets of (p, a) for both 
the prompt and afterglow cases. These are labeled as A™, 
Br9o, and C7-90 (A200, B200, and C200), and shown in Figure[2] 
In Figure |4] we show the resulting fluence distribution corre- 
sponding to each of these models. 

EGRET results imply that during the prompt emission 
phase, 0.003 < r] < 0.06. As we discussed in § 12.11 the 
low number of photons in the bursts detected by EGRET, as 
well as their spectrum, implies that the detections of prompt 
photons are most likely to have been dominated by the high- 
energy tail of the synchrotron emission; i.e., 77 k, ^ syn in Fig- 
ure |2 top). In fact, simply extrapolating synchrotron tail of 
many BATSE bursts up to ^100 MeV regime, using inferred 
values for their v sya and a-i, gives a value of 77 syn which is con- 
sistent with the one obtained here for the prompt phase. The 
harder IC prompt emission, however, can still have as much as 
10 times larger fluence than that of the synchrotron emission 
in EGRET window, without being detected. Therefore, this 
figure also sets an upper limit on the ratio of the IC and syn- 
chrotron components of 771c < 0.6, as larger 771c gives en ough 
photon fluence detectable by EGRET. As we showed in § 12.11 
theoretically we predict 771c ~ 1 .2 (for EGRET) with a canon- 
ical set of parameters. Although this appears to imply that 
the current bound from EGRET already excludes the canon- 
ical model, we cannot make such a strong statement given 
the current uncertainties of many relevant parameters. There- 
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FIG. 5. — Illustrative constraint plot on t e /eB and 7,,, from EGRET data 
for (a) prompt and (b) afterglow phases, obtained with canonical values for 
other parameters. The left and right regions of the solid curve in panel (a) is 
excluded and allowed regions respectively, and regions between two dashed 
(dotted) curves in panel (b) show allowed regions corresponding to 0.013 < 
i~i\c < 0.09 (0.006 < r)ic < 0.13). Note, however, that these regions could 
easily change depending on values of other parameters. 



fore, a more conservative statement would be that the cur- 
rent EGRET bound is barely consistent with the predictions 
of the SSC within the internal shock model. We may interpret 
the bound 771c < 0.6 as constraints on e e /es and j m , which 
is shown in Figure Ufa). As the Klein-Nishina suppression 
(£kn) becomes significant for large j m , we have only modest 
limit on e e /e/j in such a regime. However, one should keep 
in mind that these are order of magnitude constraints, which 
may farther vary with other parameters, such as v^ a ,a2 and 
IV Much better constraint plot is expected with the future 
GLAST data, where hopefully, 771c will be measured for many 



TABLE 1 
GLAST-LAT FLUENCE SENSITIVITY 



a 


to [s] 


F\im(t < to) [erg cm 2 ] 


Fj im (10 3 s) [erg cnT 2 ] 


2.3 
2.0 
1 ,0 s 


650 

650 
650 


4.5 X 10~ 7 

6.6 x 10~ 7 

5.2 x nr 6 


5.6 x 10~ 7 
8.1 x 10~ 7 
6.4 x 10~ 6 



NOTE. — Parameters of point-source flue nce sensitivity (integrated over 
30 MeV-30 GeV) of GLAST-LAT (see eq. Q_]). The power-law index is 
-a, and the unit in fluence limit F[[ m is erg cm -2 . The detection criterion for 
t < to is five photons, and significance for t > to is 5<r, where to = 650 s is the 
transition time. 

a Here we considered a detection based on the number of photons in the en- 
ergy range 30 MeV-30 GeV. A higher to and more sensitive background lim- 
ited threshold can be obtained for a = 1 if a higher energy range is considered 
(see text and Appendix IBV 

individual bursts. 

During the afterglow the synchrotron emission is much 
softer than during the prompt phase, and therefore, the IC 
component is expected to dominate EGRET observations also 
near its lower energy-band limit. Moreover, the fact that the 
number of bursts detected by EGRET during the afterglow is 
higher than the number detected during the prompt emission 
suggests that here EGRET is likely to have detected the ac- 
tual IC component of the afterglow. The spectral index of the 
GeV afterglow in EGRET window during the first 200 s is ex- 
pected to be a = 1 .5-2, implying that the evaluation of /1 in the 
bottom panel of Figure [2] which assumes a spectral index of 
-2.4, might be larger by at most a small factor (^2-3). Thus, 
for the afterglow, 771c ~ 0.01-0.1. We then compare this re- 
sult with the theoretical expectation of 771c in equation (f2]). But 
first we need to estimate the value of F syn /FM e v where F syn is 
measured during the first 200 s following Too and Fjyiev is the 
prompt emission fluence. We use the Swift GRB table 4 which 
provides X-ray afterglow fluences several tens to several hun- 
dreds of seconds after the bursts, as well as the prompt MeV 
fluences. Using only bursts where the X-ray observation starts 
after Tgo but no more than 300 s after the burst we find a dis- 
tribution of F S y n /FMev that ranges from 10~ 3 to 0.1, with the 
central value of ~10~ 2 . Thus afterglow theory with canonical 
parameters predicts 771c ~ 10~ 2 with a large spread, consistent 
with EGRET constraints. Figure |3b) shows the interpreta- 
tion of EGRET constraint on 771c (Fig.|2]i as that for e c /es and 
7 m , assuming canonical parameters and F syn /F\[ e v = 10~ 2 . Al- 
though this allowed region may change with other model pa- 
rameters, again one cannot have too large value of 7,,, because 
of the Klein-Nishina suppression factor £kn- 

4. IMPLICATION FOR GLAST 

We now move on to discussions on implications for GLAST 
using the obtained constraints on 7/ in the previous section. 
First we estimate the sensitivity of LAT on board GLAST for 
prompt and afterglow GeV emission, based on its published 
sensitivity to steady point sources, 5 which is 4 x 10~ 9 cm -2 s _1 
above 100 MeV at 5a with a power-law index of -2. This sen- 
sitivity is obtained by a one-year all-sky survey during which 
any point source is observed for ^70 d (the LAT field of view 
is 2.4 sr). 6 Therefore during the background-limited regime 
(when t is large enough that many background photons are 
observed) the flux limit scale with t as 4 x 10~ 9 cm" 2 s" 1 

4 http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table/ 

5 http://www-glast.slac.stanford.edu/ 

6 We assume here a step function for the LAT window function. 
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TABLE 2 
GRB Rate at GLAST-LAT and Contribution to the EGB Flux 



Model 


Rate at GLAST 


4gb 


[GeV 


cm ' 


s 


1 sr 1 ] 


Ar90 


15 yr 1 


6.3 x 


io-'° 


(1 + 


'/IC 


/Vsyn) 


Br90 


20 yr" 1 


8.4 x 


io- lu 


(1 + 


'/IC 


Afcyn) 


c r90 


10 yr" 1 


4.4 x 


io- lu 


(1 + 


'/IC 


/»?syn) 


A200 


20 yr" 1 




8.9 


< 10" 


111 




B200 


30 yr 1 




1.3 


x 10 


-9 




C200 


15 yr 1 




6.5 


< 10" 


10 





NOTE. — The estimate of detection rate with GLAST-LAT (for a = 2.3), 
and expected EGB intensity, for models A, B, and C of the prompt (during 
T90) and afterglow phases (during 200 s after Tgo). The correction factor 
1 + ^Tic/Vsyn for /egb in the case of prompt emission could be as large as 
~10. Also note that these estimates are quite conservative. See discussions 
in §§|4]and|5]for more details. 

(f/70 d)" 1 / 2 . During the photon-count-limited regime (when t 
is so small that less than one background photon is expected), 
in contrast, the detection limit is at a constant fluence. There- 
fore the fluence sensitivity of the GLAST-LAT detector is 



fltaCO : 



■Fiim(fo) 



1/2 



[t<t ], 
[t>t ], 



(13) 



where to = 650 s represents the time when the transition from 
photon-count-limited to background-limited regime occurs in 
the LAT case. Note that equation ( fT3l is for the limiting 
fluence, the time-integrated flux, rather than the flux. This 
limit is more natural in the photon-count-limited regime and 
it is more relevant to EGRET constraints that we derived in 
the previous section. Detailed derivation of this sensitivity is 
given in Appendix [B] In Table Q] we summarize the values 
of to and F\[ m (t) for a few cases of power law index -a and 
integration time t . The values of -fii m (f) for t <C to in the ta- 
ble are determined by criteria of five-photon detection, while 
those for t > to are by 5a significance. The fluence we argue 
here is the one integrated over 30 MeV-30 GeV, in order to 
compare with the EGRET fluence upper bounds. 

In the case of background-limited regime, it might be more 
appropriate to use higher energy threshold (instead of 30 
MeV) especially for hard source spectrum, because the back- 
ground spectrum falls steeply with frequency (a ~ 2.1). We 
may find optimal low-frequency threshold depending on spec- 
tral index of GRB emissions; it is higher for harder spectrum. 
Thus, we should be able to improve the fluence sensitivity for 
background-limited regime, compared with the figures given 
in Table Q] In addition, transition from photon-count to back- 
ground limited regime would occur later than 650 s. For our 
purpose, however, as time scales we consider (790 for prompt 
emission and 200 s after T90 for afterglows) are both during 
photon-count-limited regime, the consideration above does 
not apply and we can use full energy range (30 MeV-30 GeV 
for EGRET) to collect as many photons as possible. 

GLAST is also equipped with the GLAST Burst Monitor 
(GBM) instrument, dedicated for the detection of GRBs. It 
detects photons of 8 keV to more than 25 MeV and its field of 
view is ~8 s r. The expected rate of GRBs that trigger GBM 
is ~200 yr" 1 dMcEnery & Ritzll2006l) . which is almost as high 
as BATSE rate. Each year, about 70 out of these ~200 bursts 
should fall within the LAT field of view. Given the distri- 
bution of fluences (Fig.|4]i and the LAT sensitivity (Table [U, 
we can estimate the fraction of GRBs that would be detected 
with LAT. In Table [2] we show the expected LAT detection 
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FIG. 6. — Contribution to the EGB intensity /egb from GRBs of given 
fluence, for (a) prompt (during Tgg) and (b) afterglow emission (200 s after 
T90) phases. In each panel, three models A-C are shown. Note that for the 
prompt phase, the fluence is that for synchrotron radiation, and that for IC 
component could be even larger (see text). 



rate for a = 2.3, which is ^20 yr" 1 for the best-fit models 
of the EGRET data for both the prompt and afterglow emis- 
sions. The prompt phase estimates are for detections of the 
synchrotron component in the ^100 MeV range. Given the 
large effective area of the LAT it is expected also to detect > 
GeV photons from the IC component and identify the spec- 
tral break associated with the transition from the synchrotron 
to IC component, thereby directly testing the SSC model. 

The estimates given in Tableware fairly conservative. First, 
while we used five-photon criterion for the detection, even 
two-photon detection should be quite significant, because the 
expected background count is much smaller than one photon 
during 790 and the following 200 s that we considered. Sec- 
ond, Swift can find dimmer bursts than GBM. Although the 
discovery rate is not as high as that of GBM or BATSE, it 
would still be able to find tens of new GRBs in the LAT field 
of view. Thus the true rate would likely be larger than the 
figures given in Table [2] 

5. IMPLICATION FOR THE EXTRAGALACTIC GAMMA-RAY 
BACKGROUND 

All the GRBs except for those detected by EGRET shoul d 
contribute to the EGB flux to a certain extent (lDermerl l2007). 
This may be computed as 



7egb : 



4-TT 



dP 

dF F— [l-e(F)] : 
dF 



(14) 



where F is EGRET fluence in 30 MeV-30 GeV, dP/dF is 
the normalized distribution of EGRET fluence (eq. Ifl2ll and 
Fig. HJl, and Rgrb ~ 2 d" 1 is the occurrence rate of GRBs from 
all sky. The factor 1 - e(F) takes into account the fact that 
very bright GRBs cannot contribute to the EGB because they 
would be identified as point sources (but see discussions be- 
low). Figure|6]shows differential EGB intensity dI^c^/d\ogF 
that represents contribution from GRBs of a given fluence, 
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for prompt and afterglow phases. In the third column of Ta- 
ble 12 we show the EGB intensity due to prompt and after- 
glow phases of GRBs which is ~10~ 9 GeV cm" 2 s" 1 sr . On 
the other hand, in the same energy range, EGRE T measured 
the EGB flux to be 10" 5 GeV cm" 2 s" 1 sr" 1 (ISreekumar et al.1 
119981) . Therefore, GRBs that were detected by BATSE but 
were not detected as point sources by EGRET contribute to 
the EGB at least ^0.01%. Again, we note that the estimates 
for the prompt phase are those of synchrotron component. We 
thus need to take the predicted IC contribution into account, 
which is represented by a correction factor 1 +??ic/??syn in Ta- 
ble |2] Since this factor could be as large as ~10 according to 
the discussion in §|5] EGB flux due to prompt phase of GRBs 
could also becomes ~10 times larger, which makes GRB con- 
tribution as large as ~0. 1 % of the observations above ~GeV. 
In any case, the contributions from other astrophysical sources 
such as blazars are ex pected to be more significant than GRBs 
(e.g. lAndo et alj|2007l and references therein). 

Additional contribution to EGB is expected from a large 
number of GRBs that point away from us and therefore 
would not have been detected with BATSE. The emission 
from these bursts points toward s us once the external 
shock decelerates (Rhoads 1997). Since the total GeV 
energy emitted every decade of time during the afterglow 
is roughly constant, the contribution of these GRBs to 
EGB can be estimated by the GeV emission of the bursts 
that were detected by BATSE. Similar contribution is 
expected from bursts that points towards us but that are 
too faint to be detected by BATSE, if the GRB luminosity 
function behaves as (f>(L) oc L~ 2 a s suggested by the uni- 
versal structured jet model (Lipunov, Postnov, & Prokhorov 
20011: iRossi. Lazzati. & Reesl 120021: IZhang & Meszarosl 
2002t iP erna. Sari . & Fraill I2003L see however 

Guetta. Piran. & Waxmanl 120051) . Therefore the contri- 
bution of bursts that were not detected by BATSE to EGB 
can be estimated by the afterglow fluence of the detected 
bursts, assuming no contribution from bursts with only an 
upper limit. This is a reasonable estimate since the GeV flux 
is dominated by the few brightest bursts in GeV which are the 
most likely to be detected. Taking the fluence of the detected 
GeV bursts as the logarithmic mean of these upper and lower 
limits implies 7 E gb ~5x 10" 9 GeV cm" 2 s" 1 sr" 1 , a GRB 
contribution being ~0. 1 % of the EGB. 

Finally, we note that there is a big uncertainty in remov- 
ing the G alactic foreg round contamination from the total dif- 
fuse flux dKeshet et a l. 2004). Additionally, EGRET obser- 
vations do not constrain TeV emission that cascades down 
into the GeV range for GRBs at cos mological distances 
dCasanova et alj|2007t iMurase et alj|2007l) . Thus, if the fore- 
ground subtraction was indeed underestimated or if GRB TeV 
emission is not negligible, then GRB contribution might be 



much more significant than the estimates here. 

6. SUMMARY AND CONCLUSIONS 

The GLAST satellite would enable us to test high-energy 
emission mechanisms of GRBs. If this emission will be found 
to be consistent with SSC then its observations would con- 
strain physical parameters such as e e /e# ratio and the bulk 
Lorentz factor of the jet, IV The EGRET instrument on 
board CGRO, while less sensitive than the GLA5T-LAT de- 
tector, identified several BATSE GRBs with GeV photons. 
In addition, stringent upper limits for ^100 GRBs were put 
on fluences in the GeV b and by analyzing the EGRET data 
dGonzalez Sanchezll2005l) . 

In this paper, we further extended this EGRET result, com- 
paring with the SSC emission model. Following theoretical 
models of SSC, we assumed that there is a linear correla- 
tion between fluences in BATSE and EGRET energy bands, 
and that the proportionality coefficient r/ follows a log-normal 
distribution. We found that the predictions from the SSC 
model using canonical parameter values is fully consistent 
with EGRET fluence measurements and upper limits for both 
the prompt and afterglow phases. During the course of show- 
ing this result, we properly took the Klein-Nishina feedback 
effect into account in the theoretical calculation. The best-fit 
value of the coefficient was log 77 ~ -1 .5 for both the prompt 
and afterglow emissions, and it is already stringent enough 
to test the SSC scenario. The limits for the prompt emission 
phase are for the synchrotron radiation, and thus if we con- 
sider the IC component as well, the value of r\ could be larger 
by up to one order of magnitude. 

The obtained r\ distribution, together with the BATSE flu- 
ence distribution, gives the expected fluence distribution in 
the GeV band, which is shown in Figure |4] As the GLAST- 
LAT detector covers EGRET energy band, we can predict the 
detectable number of GRBs with GLAST from the distribu- 
tion of /"egret, given the GLA5T-LAT sensitivity. Our con- 
servative estimate using the five-photon criterion is that about 
~20 GRBs among those detected with GBM would be de- 
tected with GLA5T-LAT each year. This number could be 
even larger if we use fewer-photon criteria. The fluence dis- 
tribution can also be used to estimate the GRB contribution to 
the EGB intensity. We found that the contribution would be 
at least ~0.01% but is likely to be as large as ~0.1%. 
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APPENDIX 



KLEIN-NISHINA FEEDBACK ON HIGH-ENERGY EMISSION 

We shall find an analytic expression for £kn due to the Klein-Nishina feedback. To simplify the argument such that we can 
treat it analytically, we make the following approximations: (i) an electron with a fixed Lorentz factor 7,, radiates mono-energetic 
synchrotron photons; (ii) the same electron upscatter a given synchrotron photon to another monochromatic energy, which is 
increased by a factor of 7 2 ; (iii) vf v of both synchrotron and IC photons peaks at v syn (a synchrotron frequency corresponding to 
7 m ) and ^ic(= 7m^ S yn; if there is no Klein-Nishina suppression), respectively; (iv) the Klein-Nishina cutoff occurs quite sharply 
above its threshold; (v) both cooling and self-absorption frequencies are much smaller than the frequency region of our interest; 
and (vi) electrons cool so quickly that any dynamical effects can be neglected. With these approximations, expressions for 
the ratio of power of synchrotron and IC radiations from a given electron Y(-f e ) = P\c(je) / Psyn("fe) simplifies significantly. In 
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FIG. 7. — Ratio of IC to synchrotron power Y by an electron with Lorentz factor ■y,„ as a function of Tkn/Tm' Points represent numerical solutions of eq. IAU 
for various values of e e /es, and solid line is an analytic fit (eq. (3)). 



particular, according to the assumption (iii) above, we have (e ( ,/eB) 1 ' 2 c;KN ~ Y(^ m ). This is given as 

e e p/2-l 



T(7m) ■■ 



e B p-\ 



V^syiv 



dv —, — 



l+Y(\y>/v>-\V^ m ) 



-v +■ 



brf„. 



(Al) 



where p is electron spectral index, is the step function, and primed quantities are evaluated in the rest frame of the ejecta (e.g., 
v' = v/Ti„ where v is the frequency in an observer frame). 

A detailed derivation as well as numerical approaches are given elsewhere (Nakar, Ando, & Sari, in preparation), but at least 
this equation can be understood qualitatively. For a given electron with Lorentz factor j m , the synchrotron power does not depend 
on whether the Klein-Nishina suppression is effective or not. On the other hand, the IC power does, because it is proportional 
to the energy density of seed (synchrotron) photons integrated up to some cutoff frequency; synchrotron photons above this 
frequency c anno t be IC scattered efficiently by the electron with 7,,, because of the Klein-Nishina suppression. The integrand 
of equation (lAlt represents the synchrotron spectrum. More specifically, assuming there is no Klein-Nishina suppression, the 
spectrum is simply given by f v i oc max[z/,z; s ' ]" (p " 1) ' /2 (z/)~ 1 ' /2 ; the step function then represents the Klein-Nishina cutoff. The 
factor 1 + Y in the denominator of the integrand accounts for the suppression of the electron distribution function due to the 
enhanced IC cooling; i.e., dN e /d^f e oc (dj e /dt)~ l oc [P syn (7e) + -Pic(7e)]~ 1 oc [1 + Y(j e )]~ l . These electrons are ones that emit 
synchrotron photons of a given frequency v'. Recalling the relation % ex v' l l 2 , their Lorentz factor is given by (v'/v' n )'/ 2 7 m , 
which appears in the argument of Y in the integrand. Finally, the other constants in equation (1AU are chosen so that we have a 
proper relation for the fast cooling, Y(l+Y) = e e /e£, if w e tur n off the Klein-Nishina cutoff and have constant Y. 

Now we shall find analytic expressions of equation (IAU in asymptotic regions. We start from the case of 7,, 
m e c 2 / hv' which is equivalent to v' < m e c 2 jhr( m . The integration then becomes 



< 



7kn 



r(7m) = 



e e p/2-1 
e B p-\ 



dv' 



(<yr 1/2 



l+ytfi//^] 1 /^) 



5£ , {K yn yi 2 -\v'r pl2 



(A2) 



*V Then 



We assume that the function \+Y varies rather mildly in the integrand, so that in the argument of Y we may use v' = "' 

the integral can be evaluated analytically, and gives T(7m)[l +YQy m )] = e e /eB- When e e ^> es, we have Y(-/ m ) = (ee/tB) 1 ' 2 , which 
is the same result as in the case of no Klein-Nishina suppression. This makes sense because the condition 7,,, < 7kn indicates 
that the electrons with j m is below the Klein-Nishina threshold with seed photons at frequency v' that dominate the synchrotron 

power. On the other hand, when j m > 7kn (or v' > m e c 2 /hj m ), equation (IA1I) becomes 



Y(j m ) ■■ 



p/2-1 r^ j , (*vr 1/2 



£b p- 



1 



dv 



e e p- 



\+Y{[v'/v'Yl 2 lm ) e B p-\ 



Jm_ 

7kn 



-1/2 



1 



l+F([7m7KN] 1 / 2 )' 



(A3) 



where in the second equality, we used v' = m e c 2 /h r y m for the argument of Y. When 7„,/7kn is large enough so that 
T([7m7KN]'' 2 ) <§C 1, then equation ( IA3b immediately gives asymptotic solution for T(7,„). When j m is in the intermediate regime, 
we can still get analytic expressions, which how ever are given elsewhere because they are somewhat complicated. Here we 
simply show numerical solutions of equation (lAlb as a function of 7kn/7»i f° r various values of e e /es. We show these results as 
well as a simple fitting form (given by eq. 0) in Figure [7] Thus, equation (0 provides fairly good fit to the results of numerical 
integration of equation ( IAU . 

FLUENCE SENSITIVITY OF GLAST 

For a steady point source with a spectral index of -2, the sensitivity of GLAST-LAlT to its flux above 100 MeV is 4 x 10~ 9 
cm" 2 s _1 at 5<r significance for a one-year all-sky survey. Considering the field of view of GLA5T-LAT, 2.4 sr, this survey time 
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corresponds to 70 d exposure time to the source, and therefore, the sensitivity to the number fiuence integrated over this time 
scale is 2.4 x 10~ 2 cm" 2 . In this section, we generalize this limit to an arbitrary spectral index, -a, and exposure time, t . 

Before starting the discussion, we define the differential number and energy fluences, and integrated number and energy 
fluences (all quantities are time-integrated): 

d -^=CE-«, Z=CE l -°, (Bl) 

dE dE 



p\-a_ 

17 /--i max 



?l-a 



1- 



r?2—a rp2—a 



(B2) 



where C is a coefficient, and E m [ n and E max are the energy band boundaries. 

The fiuence sensitivity for a one year exposure is within the background-limited regime — namely within one year many back- 
ground photons are expected to be detected within the point-spread-function of the detector. In the case of GLA5T-LAT, back- 
grounds are the EGB or Galactic foreground emissions. Therefore, we start our discussion from this background-limited case. 
Let us define this background rate of GLAST 'by Nb g , for which we assume E~ 2A spectrum and use the energy-dependent angular 
resolution and on-source effective area A e ff(E). 7 The criterion of point-source detection is 



AL > ALiim = a 



v 7 



7,lim : 



where a represents significance of detection, and photon count from the source is obtained by 



al 



£n,in 



' dE^A^E). 
dE 



(B3) 



(B4) 



Therefore, using equation ( IB lb in equations ( IB4b and ( IB3I ), we can obtain the sensitivity to the coefficient Ca m as follows: 



0™=^, 



7,lim 



dEE- a A et{ (E) 



(B5) 



and then using equation (1B2| |. this can be translated into the sensitivity to the number and energy fluences, F N r lm and Fn m . 
We here note that Ci; m depends on t, a, £„,!„, and £ max , while N^n m depends only on t, E m j n , and £ m „. In this background- 
limited regi me, the time dependenc e is Fr im ex t '/ 2 from equation (IB31 I. We confirmed that, using the EGB intensity measured 
by EGRET TSree kumar et al. ( 1998) and energy-dependent angular resolution of LAT, we could obtain the limit comparable to 
F)v,iim = 2.4 x 10~ 2 cm" 2 , in the case of a = 2, t = 70 d, E^ m = 100 MeV, E m . dx = oo, and a = 5. The results of this procedure for 
several values of interest of a are summarized in equation ( fT3l and TableQ] Here, we used EGRET energy range, i.e., E^ = 30 
MeV and E m;ix = 30 GeV, but we can instead adopt different values. 

If the time scale is short such that A^ 7i ii m < 1, then the argument above does not apply, but the sensitivity is simply obtained 
by the expected photon count from the source. In this photon-count-limited regime, we can evaluate the fiuence sensitivi ty by 
requiring N 7 to be a few; here we use A^ 7 = 5. One can obtain the corresponding Cy lm by solving this criterion using equation < IB4b . 
This time, C\\ m is independent of t . Then again using equation dB2b . one can get the fiuence sensitivity in this regime as shown in 
TableQ] 
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